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Sumario. Reportamos el calculo del factor de anti-apantallamiento de Sternheimer, y, a partir de la teoria de pertur- 
I baciones de primer orden. Como base de funciones se utilizan los orbitales electronicos de Hartree-Fock, expandi- 

' O dos en estados atomicos hidrogenoideos. Las dependencias y{r), calculadas para las corazas electronicas internas de 

Ch los sistemas Fe 3+ y Cu l+ , son reportadas aqui y comparadas con valores existentes en la literatura, obtenidos por 

O metodologfas alternativas. 

Abstract. We report a calculation of the Sternheimer anti-shielding factor, y, by means of first order perturbation 
theory. In quality of basis functions, we use Hartree-Fock electronic orbitals, expanded on hydrogenic atomic states. 
The computed y(r) for Fe 3+ and Cu l+ inner electronic cores are reported and compared with literature values, obtained 
from alternative methodologies. 
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1 Introduction The next step, in the way of getting a deeper under- 

standing of the observed 57 Fe NQS data for this mate- 
In Ref. [1], the reported 51 Fe Moessbauer spectroscopy rial, must rely on an improvement of the present knowl- 
nuclear quadrupole splitting (NQS) data, measured on edge about oxygen nearest-neighborhood disorder in the 
YBa 2 CuT,0 1 - x (YBCO) samples, have been ascribed to YBCO unit-cell basal plane. This is related to the distri- 
57 Fe occupied structural sites. In this Ref., the oxygen bution of the electronic charge in the neighborhood of the 
nearest neighbour disorder, around the 57 Fe occupied Cu(l) site. 

Ck(1) crystal site in the unit cell base plane, was taken This purpose can be reached by introducing quantum 

into account. The electrical field gradients (EFG) were mechanical calculations in the framework of the Heitler 

computed according to the point charge model (PCM), and London approximation (HLA) 2 , which still offers, for 

On this ground, satisfactory qualitative assignment of all complex systems, a valid approach for EFG calculations, 

Moessbauer spectroscopy quadrupole splitting data were alternative to density functional theory (DFT). 313 

achieved for two YBCO crystalline phases: the oxygen- In the HLA, it is assumed that only valence electrons 

deficient tetragonal and oxygen-rich orthorhombic ones. are involved in chemical bonds, while the atomic inner 

electronic core (IEC) remains unperturbed, and follows 



the free-atom electronic structure. The atomic IEC po- 
larizations effects, connected to EFG calculations, are 
described through the Sternheimer anti-shielding factor 
14 ~ 17 , y(r), a kind of quadrupole polarizability coeffient. 
This y(r) can be expressed as the ratio between the IEC 
induced EFG components at the nuclear site, and the 
corresponding electrical field gradients arising from an 
external point electrical charge, placed at a distance r to 
the atomic nucleus. 

The knowledge of y(r) is essential to insure a suc- 
cessful application of the HLA for EFG calculations. 
However, there exists a wide spread of the reported y M 
(y( r )lr->°o) values 14-20 for the Fe i+ and Cu 1+ electronic 
systems. On the other hand, the full y(r) dependence have 
only been reported for the Fe inner electronic core. 19 

Therefore, for the HLA application to NQS data 
evaluation purposes, in particular to 51 Fe and 63,65 Cm hy- 
perfine interaction in the YBaiCu^Oj-^ system, further 
research work is required. The present contribution ad- 
dresses the perturbative calculation of y x and y(r) with 
the use of Hartree-Fock (HF) electronic wave functions, 
expanded on a hydrogenic basis set. 22 This methodol- 
ogy is applied to the Fe 3+ (Ne = 23) and C« 1+ (Ne = 28) 
electronic systems. 

2 Theoretical Procedure 

Hartree-Fock Calculations 

Fe 3+ and Cu 1+ inner core electronic states are de- 
fined as ls 2 2s 2 2p 6 3s 2 3p 6 3d 5 and ls 2 2s 2 2p 6 3p 6 3s 2 3d w , 
respectively . Such configurations were computed in the 
nonrelativistic HF Approximation following an algorithm 
developed earlier 22 , where each HF orbital is a linear 
combination of hydrogenic functions^: 



(i) 



where a denotes (n,l,m,m s ); a given hydrogenic quan- 
tum numbers set, and the hydrogenic wave functions \a 
are defined by: 



Xa = >f'nlmm,( r ' S -> 0- s z) = ~ P nl(r)Y bn (6 , 0)f m (s z ) , (2) 

r 

where: 



Z( B -Z-l)! /2Zr\ ,+1 f 2l+l l2Zr 
F " l[r) - in 2 [(n + iy.A n J { n 



(3) 



with respect to the quantum numbers m and m s . The re- 
sulting partial charge densities are spherically symmetric, 

Pn'l = Zm,m s If^mm, ( r > °> *> S z) I* = (t) ^if W. lead- 

ing to a vanishing contribution to the EFG tensor at the 
nucleus. 

Sternheimer anti-shielding factor calculation 

The Sternheimer anti-shielding factor y(r) is related to 
IEC polarization effects due to an external electric charge. 
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Figure 1. The unperturbed inner electronic core and the 
external point charge used in the present calculation. 

Fig.l schematically shows the inner core and the ex- 
ternal point charge q used in our calculation. The unper- 
turbed inner core charge density is spherically symmetric. 
The coordinate system origin is fixed at the atomic nu- 
cleus (nuclear dimensions are neglected), r, and R q repre- 
sent the ith inner core electron, and the external charge q 
coordinate vectors, respectively. 

The external charge induces a polarization of the elec- 
tronic cloud. The rate between Gq z , the z component of 
the EFG due to the polarized inner electronic core, and 
G qz , the z component of the EFG caused by the bare point 
charge, both taken at the nucleus site, is defined as the 
Sternheimer anti-shielding factor, according to the equa- 
tion, 



y{R q ) = 



J Sp(r, 



P2(COS0) 

0) 5 dv. (4) 



Here, 6p(r, 8) represents the polarization charge (devi- 
ation from radial symmetry). It is computed in first order 
perturbation theory from the following quadrupole pertur- 
bation Hamiltonian (in atomic units), 



£ [ -?« Ii (n 2 /Rq 3 ) Pi (cosO), for n < R q 
Q \ -qext Hi [R q 2 /n 3 ) Pi (cosO), for R q < r, 



That is: 



occupied occupied 

6p(r,ff)= J 8 Pi (r,6)= J 2^, (6) 



In the last Eq., ^^'(^f^) are Laguerre polynomials, 
Yi m (9,<p) are spherical harmonics 21 , and £ ms (s z ) are spin 
functions with defined projection along the direction OZ. 

Restricted HF calculations are performed, leading to 
quantum states, tp9, which are energetically degenerated 
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where the ipl are first-order corrections to the wave func- 
tions: 



(7) 



and the coefficients for non-degenerate levels, for exam- 
ple, are computed as By = -((p .\ff Q \<p9)/(E° - £9). Se- 
lection rules dictate that lj = Z, or /, + 2 (radial and angular 
contributions), and lead to many vanishing terms in Eq. 
(7). On the other hand, as j increases, By decreases.) 

For the Sternheimer anti-shielding factor, we get the 
following explicit expression: 

occupied 

y(R q )= £ n, (8) 

i 

where 

Radial and angular integrations in Eq. (9) can be separat- 
edly performed. The radial integrals change sign, which 
leads to alternating contributions to y,-. 

3 Results and Discussion 

The program 

A program (antish2.001), written in C++ for the GNU 
(g++) compiler, performs the Sternheimer anti-shielding 
factor estimation for a given electronic core. It uses as 
input data the restricted HF calculations for the given sys- 
tem, in particular the IEC occupied as well as the first 200 
virtual orbitals. HF orbitals are expanded in a basis of 400 
hydrogenic states. 

In our program, the radial integrals are analytically cal- 
culated. The computation of 3J-coefficients and the diag- 
onalization of the H' Q matrix is performed with the help of 
the GNU Scientific Library (GSL). 

Computational and convergence details 

Due to / quantum number selection rules, many By co- 
efficients in Eq. (7) vanish. The convergence of the series 
is relatively slow. Notice that the partial radial and angu- 
lar contributions to y{R q ) show opposite signs (see Tables 
II and IH below). 

The number of Hartree-Fock virtual orbitals is a crit- 
ical subject for convergence. In the present calculations, 
convergence is reached when using 180 or more virtual 
levels, guaranteeing relative variations lower than 0.2%, 
as can be observed in Table I for Fe 3+ y{R q = 5) calcula- 
tions. 

The computed y(R q ) values are very sensitive to small 
induced variations of the IEC wave functions. In Fig. 2, 
we show both the unperturbed and perturbed Cu l+ 3d or- 
bitals in the neighborhood of their first maximum, which 
present a relative variation of about 0.2%. This certainly 
small functional change induces a radial contribution of 
-5.301 to y(R q = 5) (see Table II below). 



Table I 


Convergence of y(R q - 5) for Fe 3+ with respect to the 


number of virtual orbitals included in the calculation. 


Number of virtual orbitals 


y(R ) 


99 


-4.84778498 


130 


-7.71152215 


140 


-7.71152215 


150 


-7.71152215 


160 


-7.71152215 


170 


-7.87306174 


180 


-7.81577288 


190 


-7.79963097 


200 


-7.79963097 



y(R q ) and joo results for Fe 3+ and Cu l+ 
Partial contributions to y{R q = 5) for Cu l+ and Fe i+ 
are presented in Tables II and III, respectively, in which, 
for comparison purposes, values reported for y^, by other 
authors are also included. The reader may find more de- 
tails in Ref. 24. 



Table II 

Partial contributions to y(R q = 5) in Cu l+ . Radial 
and angular contributions, as well as the y x values re- 
ported by other authors are included. 



Partial 


This 


Ref. 


Ref. 


Ref. 


Ref. 


Contrib. 


Work 


14 


15 


23 


18 


2p-p* 


-0.498 


-0.6 


-0.62 


-0.58 




3p-p* 


-6.580 


-6.4 


-7.9 


-13.01 




3d-d* 


-5.301 


-3.7 


-8.5 


-16.7 




Total radial 


-12.379 


-10.7 


-17.02 


-30.29 




Angular 


0.320 


2.0 


2.0 


1.19 




y{R q = 5) 


-12.059 


-8.7 


-15.02 


-29.10 


-17.37 



Table III 

The same as Tab. II for Fe +i . 



Partial 


This 


Ref. 


Ref. 


Ref. 


Ref. 


Contrib. 


Work 


17 


19 


23 


18 


2p-p* 


-0.584 


-0.724 


0.78 


-0.68 




3p-p* 


-6.463 


-8.607 


-9.98 


-5.00 




3d-d* 


-1.094 


-2.478 


-2.39 


-1.53 




Total radial 


-8.141 


-11.809 


-11.59 


-7.21 




Angular 


0.341 


0.837 


1.0 


1.04 




y{R q = 5) 


-7.800 


-10.9719 


-10.59 


-6.17 


-5.244 


(*) Transitions to excited 


orbitals with the same 


angular 





momentum. 

Notice that our results for y in Cu l+ and Fe i+ are 
within the range reported by other authors, which, on the 
other hand, show a wide spread, presumably because of 
the sensitivity of y to small variations in the polarized 
charge densities. 
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Figure 2. The IEC 3d orbital of Cu i+ in the vicinity of its 
maximum. The quadrupole field of a charge at R q = 5 a.u. from 
the nucleus causes the perturbation. 

We show in Fig. 3 normalized, y(R q )/y m results for 
the studied Cu l+ and Fe 3+ systems. Notice that, in both 
cases, y(R q > 5) ~ y^. The dotted lines in Fig. 3, rep- 
resent fitting curves according to y(R q )/yoo = 1 + A. * 
expi-Rg/p), where A(Cu 1+ ) = -3.797, p(Cu l+ ) = 0.756, 
and A(Fe 3+ ) = -7.027, p(Fe 3+ ) = 0.550. 
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Figure 3. Cu [+ and Fe 3+ Sternheimer anti-shielding factor 
dependence on R q . Dotted lines are fitting curves according to 
y(R q )IVoo = \ + A* e { - R " lp \ where jUCu l+ ) and y m {Fe +3 ) are 
taken from Tables II and III, respectively. 

4 Conclusions 



The Sternheimer anti-shielding factor, y(R q ), was suc- 
cessfully computed in the framework of first order per- 
turbation theory and a hydrogenic function basis set. The 
HF calculations were performed in a restricted spherically 
symmetric approach. 

The computed y(R q ) dependences and y M values for 
the Cu l+ and Fe 3+ are within reasonable ranges. The cal- 
culation methodology developed in the present paper may 
be extended and applied to other atomic systems. 

The obtained y(R q ) may be used in the computation of 
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nuclear quadrupole splittings in the solid state systems. 
This is left for a future work. 
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